Steady-State Source Flow in Heterogeneous Porous Media

The flow of fluids in heterogeneous porous media is modelled by regarding the hydraulic conductivity as a stationary random space function. The flow variables, the pressure head and velocity field are random functions as well and we are interested primarily in calculating their mean values. The latter had been intensively studied in the past for flows uniform in the average. It has been shown that the average Darcy's law, which relates the mean pressure head gradient to the mean velocity, is given by a local linear relationship. As a result, the mean head and velocity satisfy the local flow equations in a fictitious homogeneous medium of effective conductivity. However, recent analysis has shown that for nonuniform flows the effective Darcy's law is determined by a nonlocal relationship of a convolution type. Hence, the average flow equations for the mean head are expressed as a linear integro-differential operator. Due to the linearity of the problem, it is useful to derive the mean head distribution for a flow by a source of unit discharge. This distribution represents a fundamental solution of the average flow equations and is called the mean Green function G _d (x). The mean head G _d(x) is derived here at first order in the logconductivity variance for an arbitrary correlation function (x) and for any dimensionality d of the flow. It is obtained as a product of the solution G _d ⁽⁰⁾(x) for source flow in unbounded domain of the mean conductivity K _A and the correction _d (x) which depends on the medium heterogeneous structure. The correction _d is evaluated for a few cases of interest.Simple one-quadrature expressions of _d are derived for isotropic two- and three-dimensional media. The quadratures can be calculated analytically after specifying (x) and closed form expressions are derived for exponential and Gaussian correlations. The flow toward a source in a three-dimensional heterogeneous medium of axisymmetric anisotropy is studied in detail by deriving ₃ as function of the distance from the source x and of the azimuthal angle . Its dependence on x, on the particular (x) and on the anisotropy ratio is illustrated in the plane of isotropy (=0) and along the anisotropy axis ( = /2).The head factor k ^* is defined as a ratio of the head in the homogeneous medium to the mean head, k ^*=G _d ⁽⁰⁾/G _d= _d ^–1. It is shown that for isotropic conductivity and for any dimensionality of the flow the medium behaves as a one-dimensional and as an effective one close and far from the source, respectively, that is, lim _x0 k ^*(x) = K _H/K _A and lim _x k ^*(x) = K ^efu/K _A, where K _A and K _H are the arithmetic and harmonic conductivity means and K ^efu is the effective conductivity for uniform flow. For axisymmetric heterogeneity the far-distance limit depends on the direction. Thus, in the coordinate system of (x) principal directions the limit values of k ^* are obtained as . These values differ from the corresponding components of the effective conductivities tensor for uniform flow for = 0 and /2, respectively. The results of the study are applied to solving the problem of the dipole well flow. The dependence of the mean head drop between the injection and production chambers on the anisotropy of the conductivity and the distance between the chambers is analyzed.     

杭州大剧院风压分布的风洞试验研究

杭州大剧院外观造型新颖独特，由一个近似为椭球的双曲面屋盖和一个倒置的圆锥面玻璃幕墙相交形成大剧院的主体建筑的外轮廓．本文结合该主体结构的抗风设计要求，制作了1：200刚性模型进行结构表面脉动风压平均风压分布的风洞试验研究，详细介绍了风洞试验及数据处理方法，讨论了在不同风向角下椭球面屋盖和圆锥面玻璃幕墙的风压分布．为了便于整体结构设计，本文根据结构的实际情况对椭球面和圆锥面进行分块并给出相应的分块风压数据，最后根据各风向角下的风压分布特征给出了整体结构抗风设计的最不利风向角．     

光纤陀螺捷联系统标定测试软件的设计与开发

基于 windows 平台,采用 VC 6.0 开发了光纤陀螺捷联惯组标定测试软件,着重阐述了该软件的设计思想和关键技术的实现。该软件实现了对光纤陀螺捷联组合误差系数的标定与测试,为光纤陀螺捷联系统标定测试提供了方便,具有很强的工程实用价值。     

Direct Numerical Simulation of Turbulence in a Stably Stratified Fluid and Wave-Shear Interaction

Turbulence decay in a strongly stratified medium is simulated by a direct pseudo-spectral code solving the three-dimensional equations of motion under the Boussinesq approximation. The results are compared to non-stratified simulations results. We focus on the production of mean shear energy observed in the stratified case. We then simulate the decay of stratified turbulence when affected by an initial horizontal mean flow and show that this mean flow is the major component remaining at large t. Next, we give some analytical elements on wave-shear interaction by using a simple refraction calculation with WKB hypothesis. This calculation is illustrated by simulating the interaction between one monochromatic internal wave and a vertical shear profile. We conclude that the existence of singularities in the mean shear production term in the presence of internal gravity waves may be one of the possible mechanisms involved within stratified turbulent shear flows.     

三轴转台轴系相交度误差测量方法的研究

—本文提出了一种新的轴系相交度误差测试方法。该方法不存在测量基准安装偏心问题，更适合于大型转台轴系相交度误差的测量     

船用惯性导航系统姿态测量误差辨识及其补偿方法研究

本文分析了某型号船用惯性导航系统（ＳＩＮＳ）的试验数据，结果表明，ＳＩＮＳ的姿态误差主要来源于水平测姿传感器。据此，本文对该ＳＩＮＳ水平姿态测量误差的辨识及其补偿方法进行了研究，并建立了水平测姿误差的回归数学模型，以提高ＳＩＮＳ姿态精度     

CONVERGENCE OF SIMPLIFIED AND STABILIZED MIXED ELEMENT FORMATS BASED ON BUBBLE FUNCTION FOR THE STOKES PROBLEM

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｏｒｄｅｒｔｏｇｅｔｒｉｄｏｆｔｈｅｒｅｓｔｒｉｃｔｉｏｎｔｏｄｉｓｃｒｅｔｅＢａｂｕｓｋａ＿Ｂｒｅｚｚｉ’ｓｃｏｎｄｉｔｉｏｎｓｉｎｒｅｓｅａｒｃｈｅｓｏｆｍｉｘｅｄｆｉｎｉｔｅｅｌｅｍｅｎｔ (ＭＦＥ)ｍｅｔｈｏｄｓ,ｏｎｅｕｓｕａｌｌｙｕｓｅｓｔｈｅｓｔａｂｉｌｉｚｅｄＭＦＥｍｅｔｈｏｄｓｆｏｒｍａｔ[1] .Ｔｈｅｒｅａｒｅａｎｕｍｂｅｒｏｆｓｔａｂｉｌｉｚｅｄｍｅｔｈｏｄｓ (ｓｅｅ ,ｅ.ｇ .,Ｒｅｆｓ.[1- 7]ａｎｄｔｈｅｒｅｆｅｒｅｎｃｅｓｔｈｅｒｅｉｎ) ,ｗｈｉｃｈｈ…     

ADAPTIVE MESHLESS METHOD BASED ON LOCAL FIT TECHNOLOGY

An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement.New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.     

中心裂纹圆盘应力强度因子的测试误差分析

本文在中心裂纹圆盘应力强度因子解析解的基础上,利用一阶微分法则,给出了与裂纹相对长度和加载角相关的应力强度因子(K 和K )的4个误差传递函数。这4个误差传递函数关于裂纹相对长度和加载角均是非线性的,它们既是误差分析的基础,又是合理确定裂纹相对长度和加载角的基础。分析结果表明,加载角的误差Δθ除了对纯 型K 的误差几乎没有影响,对纯 型K 影响较小外,对复合型K 、K 的误差均有较大影响。最后,本文建议裂纹相对长度的取值范围为0.4～0.6;还建议在复合型断裂试验时,必须依据对K 、K 的总体精度要求来严格控制加载角的精度。     

Multigrid Models of Composite Materials of Irregular Structure with a Small Filling Ratio

This paper considers composites consisting of a set of typical composite multigrid finite elements whose structures are regular and different. Mean local errors are proposed for multigrid modeling of composites.